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Dynamics of brane- world models of dark energy is reviewed. We demonstrate that simple dark 
energy brane models can be represented as 2-dimensional dynamical systems of a Newtonian type. 
Hence a fictitious particle moving in a potential well characterizes the model. We investigate the 
dynamics of the brane models using methods of dynamical systems. The simple brane- world models 
can be successfully unified within a single scheme - an ensemble of brane dark energy models. 
We characterize generic models of this ensemble as well as exceptional ones using the notion of 
structural stability (instability). Then due to the Peixoto theorem we can characterize the class of 
generic brane models. We show that global dynamics of the generic brane models of dark energy is 
topologically equivalent to the concordance ACDM model. We also demonstrate that the bouncing 
models or models in which acceleration of the universe is only transient phenomenon are non-generic 
(or exceptional cases) in the ensemble. We argue that the adequate brane model of dark energy 
should be a generic case in the ensemble of FRW dynamical systems on the plane. 

PACS numbers: 98.80.Bp, 98.80.Cq, ll.25.-w 



I. INTRODUCTION 



Recent observations of distant supernovae type la (SNIa) [l|, [2[ as well as other current observations of cosmic 
microwave background radiation (CMBR) anisotropics indicate that our universe is almost flat and stay in the 
accelerated phase of expansion [3|, l4[. Principally there are two alternatives of explanation properties of current 
universe. In the first approach the universe is accelerating due to presence of mysterious energy called dark energy 
X, of unknown origin, which violates the strong energy condition px +3px > 0, where px and px are energy density 
and pressure of dark energy, respectively. Alternative to the dark energy idea is the explanation of acceleration of the 
universe based on some modification to the Friedmann-Robertson- Walker (FRW) dynamics arising from new physics 
(i.e., extra dimensions). In this approach dark energy is the manifestation of modified gravitational dynamics of 
the 4-dimensional brane leading to a self-accelerated universe. In this conception which is called brane dark energy 
new physics mimic on the phenomenological level dark energy or effects of extra dimensions manifest themselves as 
a modification to the Friedmann equation which govern the evolution of the 4-dimensional brane localized in the 
higher-dimensional bulk space. In contrast to the classical Kaluza-Klein theories (when extra dimensions, which 
manifest themselves at high energies, are either compact and have finite volume), theories with infinite volume of 
extra dimensions modify the late time cosmological evolution of the universe and gravitational dynamics is preserved 
at short distances. 

The hierarchy problem [H, @, 0] and the problem of present acceleration of the Universe [1, 0] also motivate all 
theories with large extra dimensions. In the literature of the subject there are many different propositions of dark 
energy of brane origin (for review see 10]). Different brane- world scenarios were tested by recent astronomical 
measurements of distant supernovae la [11] as well as other observations from Type lib radio galaxies and X-ray gas 
mass fraction [l2|, fl3|. T he stringent constraint on the model parameters can be also obtained from baryon oscillation 
peak measurements [14], [l5|, [l6| . 

The main aim of this paper is to investigate dynamics of the most popular brane dark energy models using methods 
of dynamical systems. We investigate a class of the simplest brane models in which our observable universe is a 
(3 + l)-dimensional brane and a 3-space is homogeneous and isotropic. The physical space is embedded in a (4 + 1)- 
dimensional space called a bulk space. For all these models we obtain an analog of the classical Friedmann equation 
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with some additional terms which are functions of both the scale factor a and energy density of matter on the brane 
(see Table |Tj). If we assume that the energy density satisfies the conservation condition p = —3H(p + p), where p — wp 
is the equation of state, then p = p(a) and the modified FRW equation on the brane can be represented by 



where H = (In a)' is Hubble's function, k is the curvature index, p e s = p e ff(&) is the effective energy density. For the 
concordance ACDM model p e ^ = Pm,o(<V&o) -3 + A, p e ^ = — A. Without lost of degree of generality one can choose 
ao = 1 - as the present value of the scale factor; 1 + z = do/ a, where z is redshift. 

Equation (pQ) can be rewritten to the new form analogous to the first integral of the energy for a unit mass particle 
moving in the potential well 

y + V» = 0, (2) 

where a dot denotes the differentiation with respect to the rescaled time t — > r : \Ho\dt = dr and V(a) = — ^{^ m ,o& -1 + 
Qk,o + ^A,o^ 2 }; ^m,o = ff^r, ^A,o = 3^2 are density parameters for matter and the cosmological constant on the 

brane, V(a) = -^f^. 

Equation ([2j) reduces the problem of dynamics of a cosmological model to the standard problem of classical mechanics 
- motion of a particle of a unit mass in the potential well V(a) on the distinguished "zero energy level". 

For the FRW brane cosmological models there is a counterpart of the Friedmann first integral ([2]) in the form 

d 2 

— + V(a) + ^brane(a) = 0, (3) 



where V hrane (a) = -fa^ y( a ) 



6 ' r ^> 6 " 

Both relations ([2]) and (|3J) play the role of a first integral of motion for the system 



X Z\_v (4) 

y ~ dx-> 

which has the form of a 2-dimensional dynamical system of a Newtonian type and x = a/a$. 
The classical motion of the system is restricted to the domain admissible for motion: 

V = {x: V<0}. (5) 

The evolution of the system has been visualized in 2-dimensional phase space (x, x). There are two types of solutions 
of (@J): 1) a singular one which corresponds to the critical points of system (j4]) such that yo = 0, and (dV/dx) Xo = 0; 
2) a nonsingular one which lies on the algebraic curves determined by the Friedmann first integral y 2 /2 + V(x) = 0. It 
is convenient to extract the curvature term from the potential function. Then we can regard evolution of the system 
on a constant energy level E = ^Qk,o- 

An advantage of dynamical system methods is a investigation of all admissible solutions for all initial conditions. A 
picture of all evolutional paths represented by phase curves and critical points form phase portrait - global visualization 
of dynamics. The trajectory of the flat model Qk,o = divides phase portrait in to two subsets occupied by open 
models (Clk,o > 0) an d by closed models (Clk t o < 0). 

The character (type) of the critical points is, following the Hart man- Grobman theorem, determined from the 
linearized system 

(x - x )' = y, 2 ^ 

around the critical point (which is a static critical point in any case in a finite domain of the phase space). Therefore 
the characteristic equation of the linearization matrix is 

A 2 + detA = 0, (7) 

where tr A = 0, det A = V xx (xq). 

From equation ([7|) we obtain that only two types of critical points are admissible for the system of a Newtonian 
type, namely saddles if V xx (xq) < or centers in opposite case of V xx (xq) > 0. The case of V xx (xq) = (a inflection 



TABLE I: Classification of simple dark energy brane models in terms of H(x) relation and Vb r; 



-O^f^ = -in brane (a)a 2 . 
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Dvali- Turner model [171 ] 

Deffayet-Dvali-Gabadadze 
model p — wp [18| | 

Shtanov-Sahni model w — 
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point) is a degenerated case. In the first case eigenvalues of the linearization matrix are real of opposite sings. In the 
case of centers (non-hyperbolic critical points) eigenvalues are purely imaginary. 

For the conservative system it is useful to develop methods of qualitative investigations of differential equations 
[2ll | . The main aim of this approach is to construct phase portraits of the system which contain global information 
about the dynamics. The phase space (x, y) offers a possibility of natural geometrization of the dynamical behavior. 
It is in a simple 2-dimensional case structuralized by critical points or non-point closed trajectories (limit cycles) 
and trajectories joining them. Two phase portraits are equivalent modulo homeomorphism preserving orientation of 
the phase curves (or phase trajectories). From the physical point of view critical points (and limit cycles) represent 
asymptotic states (equilibria) of the system. Equivalently, one can look at the phase flow as a vector field 



whose integral curves are the phase curves. 

Thanks to Andronov and Pontryagin [22[ , the important idea of structural stability was introduced into the ensem- 
ble [23[ of all dynamical systems. A vector field, say /, is a structurally stable vector field if there is an e > such that 
for all other vector fields g, which are close to / (in some metric sense) ||/ — g\\ < e, / and g are topologically equiv- 
alent. The notion of structural stability is mathematical formalization of intuition that physically realistic models in 
applications should posse some kind of stability, therefore small changes of the r.h.s. of the system (i.e. vector field) 
doesn't disturb the phase portrait. For example motion of pendulum is structurally unstable because small changes of 
vector field (constructed from r.h.s. of the system) of a friction type g = [y, —dV/dx + ky] T dramatically changes the 
structure of phase curves. While for the pendulum, the phase curves are closed trajectories around the center, in the 
case of pendulum with friction (constant) they are open spirals converging at the equilibrium after infinite time. We 
claim that the pendulum without friction is structurally unstable. Many dynamicists believe that realistic models of 
physical processes should be typical (generic) because we always try to convey the features of typical garden variety 
of the dynamical system. The exceptional (non-generic) cases are treated in principle as less important because they 
interrupt discussion and do not arise very often in applications [24[ . 

In the 2-dimensional case, the famous Peixoto theorem gives the characterization of the structurally stable vector 
field on a compact two dimensional manifold [25[ . They are generic and form open and dense subsets in the ensemble 
of all dynamical systems on the plane. If a vector field / is not structurally stable it belongs to the bifurcation set. 

The space of all conservative dark energy models can be equipped with the structure of the Banach space with the 
C 1 metric. 

Let V\ and V2 be two dark energy models. Then C 1 distance between them in the ensemble is 



where E is a closed subset of configuration space. Of course, ensemble of all dynamical systems of Newtonian type 
on the plane is infinite dimensional functional space and the introduced metric is a so-called Sobolev metric. 

While there is no counterpart of the Peixoto theorem in higher dimensions it can be easy to test whether planar 
polynomial systems, like in considered case, have structurally stable phase portraits. For this aim the analysis 
of behavior of trajectories at infinity should be performed. One can simply do that using the tools of Poincare 
S 2 construction, namely by projection trajectories from the center of the unit sphere S 2 = {(X,Y,Z) G R 3 : 
X 2 + Y 2 + Z 2 = l} onto the (x, y) plane tangent to S 2 at either the north or south pole. 

The vector field / is structurally unstable if: 

1. there are non- hyperbolic critical points on the phase portrait, 

2. there is a trajectory connecting saddles on the equator of S 2 . 

In opposite cases if additionally the number of critical point and limit cycles is finite, / is structurally stable on S 2 . 

Let us consider 2-dimensional dynamical system of a Newtonian type ([5]). There are three cases of behavior of the 
system admissible in the neighborhood of the critical point (xo, 0) : — dV/dx\ XQ = 0: 

• If (#o,0) is a strict local maximum of V(x), it is a saddle point; 

• If (#o, 0) is a strict local minimum of V(x), it is a center; 




(9) 



• If (#o,0) is a horizontal inflection point of V(x), it is a cusp. 
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FIG. 1: Three possible types of behavior in the neighborhood of a critical point (from left to right): a saddle, a center and a 
cusp. 



All these cases are illustrated in Fig. [TJ 

It is a simple consequence of the fact that characteristic equation for linearization matrix at the critical point 
(tr A = 0) is A 2 + detA = 0, where det A = d 2 V/dx 2 \ Xo . Therefore the eigenvalues are real of opposite sign for saddle 
point, and for centers - if they are non- hyperbolic critical points - purely imaginary and conjugated. 

In Fig. [2] the phase portrait for the ACDM model is presented on compactified projective plane by circle at infinity. 
Of course it is structurally stable. Therefore, following the Peixoto theorem, it is generic in the ensemble M. of all 
dark energy models with 2-dimensional phase space because such systems form open and dense subsets. 

It is also interesting that the phase space of dark energy model can be reconstructed from SN la data set and is 
topologically equivalent to the ACDM model. Fig. [3] represents the potential function 



V(a(z)) = --(l + z) 2 



d d L (z) 
dz 1 + z 



(10) 



reconstructed from the relation (Il(z) - the luminosity distance 4 as a function of redshift z: 1 + z 
reconstruction is possible due to the existence of the universal formula 



d L (z) 
1 + z 



dz' 



Such a 



(ii) 



for the flat model. 

The main aims of presented discussion is searching for generic cases of brane- world models which are generic in the 
ensemble of brane- world models. We argue that adequate brane models should be located in the near ^-neighborhood 
of the ACDM model which is structurally stable. Therefore the global structure of their dynamics equivalent to the 
ACDM model is required. This implies that structural stability becomes requirement which contain model parameter. 



II. BRANE DARK ENERGY MODELS AS A DYNAMICAL SYSTEMS. CONCLUDING REMARKS 

Different brane models which offer explanations of acceleration of the current Universe can be characterized in 
terms of an additional term Vbrane(^) as a function of the scale factor. In Table [J we complete most popular brane 
models. They might provide answers to a number of cosmological puzzles including the issue of dark energy or type of 
the cosmological singularity. Model called Sahni-Shtanov was introduced in [20] . There are two classes of such models 



FIG. 2: The phase portrait for the ACDM model. In the phase portrait we have a single saddle point in the finite domain and 
four critical points located on the circle at infinity. They represent an initial singularity (x = 0,x = oo) or a de Sitter universe 
(deS). The trajectory of the flat model (k = 0) divides all models in to two disjoint classes: closed and open. The trajectories 
situated in region II confined by the upper branch of the k = trajectory and by the separatrix going to the stable de Sitter 
node deS+ and by the separatrix going from the initial singularity to the saddle point S correspond to the closed expanding 
universe. The trajectories located in regions I and II (which corresponds to the open universes) are called inflectional. Quite 
similarly, the trajectories situated in region III correspond to the closed universes contracting from the unstable de Sitter 
node towards the stable de Sitter node. The trajectories running in this region describe the closed bouncing universes. The 
trajectories located in the region between x axis and separatices of the critical point S correspond to the oscillating closed 
universes expanding from the initial singularity located at (x = 0, x = oo) towards the final singularity at (x = 0, x = — oo). 



Brane I and Brane II. Model called Shtanov-Sahni which describe bouncing model with timelike extra dimensions was 
introduced in [l^ . 

The Einstein equations on the brane, in general, form a very complicated system of non-linear partial differential 
equations. However, the majority of most interesting models from the cosmological point of view belong to the class 
of homogeneous and isotropic ones, for which this complicated system of equations reduces to the form of a dynamical 
system. Hence in investigation of dynamics of these models the methods of dynamical systems seems to be natural. 
A main advantage of this methods is that we can obtain on the phase portraits global visualization of the dynamics, 
i.e., all evolutional paths admissible for all initial conditions. Hence we obtain asymptotic states of the system and 
we can analyze their type of stability (character of a critical point). The full knowledge of dynamics require analysis 
of dynamical behavior not only in a finite domain of the phase space but also at infinity. In our paper such an 
analysis was performed using the technique of Poincare compactification of a plane by a circle at infinity. Although 
the application of qualitative methods in investigation of differential equations does not reveal unexpected properties 
of the brane models unnoticed until now, this methods offers a possibility of their classification in the terms of the 
evolutional paths. Moreover one detect in the phase space all global attr actors and their insets. From the physical 
point of view it is interesting to know how large is inset or outset of limit states because the probability of an initial 
state of the observation to evolve asymptotically to a limit set is proportional to the volume of its inset. For example 
attractors have open insets, so they are probable (the inset of an attractor is called its basin). 

In our analysis we mainly concentrate on the analysis of brane dynamics and its reduction to the dynamical system 
of a Newtonian type. But we also note that some alternative representation of the dynamics can be useful in the case 
of Dvali and Turner model. In some sense in this paper we extend previously introduced methods [U, [28[ in large 
class of brane models. 

Apart of detailed analysis of the brane dynamics using dynamical systems techniques we are looking for the generic 
dynamical systems of brane origin. Usually exceptional cases are more complicated from the mathematical point of 
view and they interrupt the discussion. Moreover dynamicists believe that they should not arise very often in the 
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FIG. 3: The potential function as a function of the scale factor expressed in its present value ao = 1 for the reconstructed best 
fit model is given by the solid line. The confidence regions la and 2a are drawn around it. The phase portrait obtained from 
this potential (best fit) is equivalent to ACDM (see Fig. [2]). The value of the redshift transition estimated from SNIa data 
(Gold sample) is about 0.38 (see |26j]V 




FIG. 4: Structurally stable phase portraits of Deffayet-Dvali-Gabadadze model model with ^ m ,o = 0.3, Q rc = 0.15 and: a) 
w = 0, b) w = 1/3. Dotted trajectories represent flat model trajectories Qk,o = 0. 
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a) 




FIG. 5: Degenerated phase portraits of Deffayet-Dvali-Gabadadze model model with n m ,o 
b) w = — 1 and c) w = —4/3. 



0.3, Q rc = 0.15 and: a) w = -1/3, 



application, because they are exceptional (|24|, p.349]). Abraham and Shaw pointed out in their famous book that 
a considerable portion of the history of the mathematical dynamics has been dominated by the search for generic 
properties. The idea was to characterize a class of phase portraits that are far simpler than arbitrary ones. The main 
motivation of the search is to complete classification of the complexity of the phase portraits (modulo exceptional 
cases). 

This idea was achieved for dynamical systems on the plane by Peixoto due to fundamental results obtained by 
Russian mathematicians Andronov and Leontovich. Peixoto theorem characterize generic dynamical systems on the 
plane in tools of notion structural stability which was previously introduced by Andronov and Leontovich. This 
theorem states that structurally stable systems on the 2-dimensional closed space forms open and dense subsets in the 
functional space of the dynamical systems on the plane. The structurally unstable dynamical system are exceptional 
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FIG. 6: Sahni-Shtanov model for ^ m ,o = 0.3, ^a,o = 0.05, Qz 5 o = 0.2 and: a) brane I model (+), b) brane II model (-). 
in this space. 

The main goal of the search of structural stability of dynamical systems of brane origin on the compactified Poincare 
sphere is to narrow down the complexity of the portraits enough to allow classification of generic cases. This was 
simply achieved because for dynamical systems on the plane with polynomial right hand sides we have simple test 
of structural stability (instability). The notion of structural stability is intuitively very simple: A vector field has 
the property of structural stability if all delta perturbations of it (sufficiently small) have epsilon equivalent phase 
portraits. The equivalence of the phase portraits is established by homeomorphism preserving the orientation of the 
phase curves (topological equivalence). 

It is pointed out that the brane world models have certain features that distinguish them form the other models of 
dark energy. They can allow for a "transient" universe acceleration phenomena which is preceded and also followed 
by matter domination era. Also brane models admit strange type of singularities which are known as a "quiescent" 
cosmological singularities. It is also possible that initial singularity can be replaced by characteristic bounce. The 
question is: Whether these phenomenas are attributes of a generic brane dynamical systems? Our answer is No. 
Typical brane models (on the plane) have global dynamics represented by phase portraits which should be topologically 
equivalent to the ACDM model. 

The brane models which a topological structure of the phase space is equivalent to the ACDM phenomenological 
model provides a simple alternative to explanation of accelerating expansion of the current universe. In the concor- 
dance ACDM model parameter A is constant with respect to redshift z and matter density evolves with redshift like 
(1 + z) 3 . Hence appear basic problem: Why do they approximately equal each other now ? This problem called 
coincidence conundrum is not solved in the frame of the A model. Alternatively there is a possibility that there is 
no dark energy, but instead an infrared modification of the Friedmann equation on very large scales. To solve this 
problem in the frame of the Deffayet-Dvali-Gabadadze (DDG) model it is required r c ~ Hq . 

We show that the DDG(+) (Fig. [4]) model has phase portrait equivalent to the ACDM model which is a purely 
phenomenological theory offering the description of acceleration of the current universe rather than its understanding 
(note that they does not explain why the vacuum energy does not gravitate |29[). 

The notion of structural stability allows us to formulate some general assumptions about any adequate brane models 
offering the explanation of acceleration of the current universe in terms of the potential function: 

• the potential function has a maximum at some finite value of the scale factor or redshift; it is a redshift transition 
value corresponding to switching a deceleration phase to an accelerating one, 

• at the late time the potential function behaves like V oc —a 2 which guarantee the existence of deS + global 
attractor. 
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FIG. 7: Shtanov-Sahni model for Q m ,o — 0.3, Qa,o — 0.45, Qdr,o — 0.15 and: a) w = 0, e = 1, b) w = 0, e = —1, c) k; = 1/3, 
e = 1, d) w = 1/3, e = —1. 

In our approach the dynamics is reduced to the 2-dimensional phase plane compactified by a circle at infinity. If 
the set of all vector fields / G C r (M 2 ) (r > 1) having a certain property contains an open dense subsets of C r (M 2 ) 
then the property is called generic. 

In the zoo of brane-world models all discussed previously phenomenological models can be recovered. In this way 
we can find counterparts of phantom models within class of Sahni-Shtanov models as well as the Cardassian one 
within the DDG(-) (Fig. [5]) brane models. 

Most interesting brane-world models should be, similar to the ACDM model, structurally stable. Such a postulate 
is satisfied by the DDG(+) (Fig. [4j) models as well as some Sahni-Shtanov models (Fig. [7j) and Dvali- Turner models 
(Fig. [8]). In our opinion they should be treated seriously as candidates for explanation of current SNIa data and other 
concordance astronomical observations. The advantage of this models is that they offer some physical mechanism of 
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b) 




FIG. 8: Phase portraits for Dvali- Turner model for ^ m ,o = 0.45. From Table |T] for a = 1 we receive a quadratic equation for 
the variable y — (H/Ho)(a/ao) and hence two solutions for the potential function: a) with a plus sign, b) with a minus sign. 




FIG. 9: a) Phase portrait for Sahni-Shtanov brane I model with £l m ,o = 0.2, Qa,o = and f^o = 0.32. Dashed lines denote 
flat model trajectories Qk,o = 0. Solid line (blue in el. version) denotes maximal allowed value of scale factor a for which 
square root in case IV (see Table UJ) is positive [27[ . b) Phase portrait for Shtanov-Sahni model with Q m ,o = 0.3, ^dr,o = 0.47, 
Qa,o = 0.12 and w = 0, e — —1. This class of brane-world models allows for a transient acceleration of the universe which is 
proceeded and followed by a matter domination era. Such models admit quiescent singularities. Model with time-like extra 
dimension can also avoid cosmological singularity by a bounce. 
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this acceleration in contrast to the ACDM model. 

From the physical point of view it is interesting to know whether certain subsets of ensemble contain an open and 
dense subsets because it means that this property is typical. 

Recently the Bayesian Information Criterion (BIC) was applied in the context of finding (choosing) adequate accel- 
erating model of the universe [30[. The authors showed preferences for the models beyond standard FRW cosmology 
whose best fit parameters reduce them to the cosmological constant model. We characterize exotic cosmological 
models in the terms of structural stability and equivalence of phase portraits. 

It was pointed out [3l| that concepts of rigidity and fragility formulated in terms of structural stability may be 
very useful for cosmology which operates by models. Moreover, as put forward many years ago by Andronov and 
Pontryagin [22[ the structural stability seams to be important condition for building adequate models of physical 
phenomenons in all branches of science. In general, there is a widespread opinion that physically realistic models 
of the world should possess some kind of structural stability. Because, there can by many dramatically different 
mathematical models all agreeing with observations and it would be fatal for the empirical methods of science [32[ 
if we are taking into account final errors of measurements [33]. In our opinion structural stability of dark energy 
models (note that ACDM is favored by astronomical data by Bayesian selection framework [34], [35[ [H, [37|) reflects 
their flexibility with respect the data fitting. 

We always try to convey the features of typical garden variety of dynamical system accelerating cosmological models. 
The exceptional cases are usually more complicated and they in principle interrupt the discussion. This prejudice 
is shared by all dynamicists ([24], p. 149]). We find that brane models in which instead of initial singularity there is 
present bounce are structurally unstable because of the center on the phase portraits. Therefore the bounce is not 
generic property of the evolutional scenario of accelerating brane models. Such a models require presence of double 
accelerating phases during the cosmological evolution. In terms of potential function of dynamical system it means 
that the existence of minimum before the maximum. If we consider the Sahni-Shtanov brane II model then for some 
set of parameters appears a quiescent singularity at which the density, pressure and the Hubble function remain finite, 
while all the Riemannian invariants diverge to infinity. Similar models belong to the class of non-generic cosmological 
models. Note that all models within the accelerating phase, which is only a transitional phenomenon, are non-generic 
too. As it was demonstrated by Sahni and Shtanov [38] there is a class of loitering brane- world models. In the phase 
portraits the corresponding evolutional paths are situated near the separatrices of saddle point (inflectional models). 
This property is generic in contrast to a quiescent future singularity. 
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